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1.  Introduction 

The  dynaaic  environments  to  which  weapon-vehicle  systems,  such 
as  those  modeled  by  Andrews  Cl],  are  subjected  include  steady  state 
harmonic  excitation,  shock  and  random  excitation.  These  occur  as  a 
result  of  external  stimuli  such  as  atmospherically  borne  disturbances 
(wind,  wake  acoustic  noise,  rotor  tip  vortex  loading)  and  through 
various  weapon-vehicle  interactions.  In  the  problems  of  this  project 
random  excitation  is  the  rale  rather  than  the  exception.  Thus  this 
report  susmarises  the  probability  techniques  necessary  for  and  their 
application  to  the  development  of  analytic  methods  for  obtaining  the 
response  of  linear  elastic  structures  to  certain  claeses  of  random 
excitation. 

This  random  vibration  response  analysis  employs  the  normal 
modes  of  a  lumped  parameter  representation  of  a  complex  system.  The 
random  forcing  functions  at  each  nods  (in  terms  of  expected  value 
and  power  spectral  density)  at  a  transformed  to  a  sat  of  modal  forcing 
functions.  Than  the  response  of  each  mods  to  a  random  forcing  function 
can  be  obtained  using  tbs  nodal  transfer  functions*  Finally  ‘•he  aodal 
responses  are  transformed  back  to  tha  physical  plana.  The  results  are 
the  statistical  expected  values  (mean,  root  mean  square,  power  spectral 
density)  of  tha  displacements,  velocities  and  accelerations  of  tha 
physical  system. 


2.  Basic  Probability  Concepts 
(a)  Random  Variable 


Let  x  be  a  real  number  resulting  from  a  measurement  and  define 
an  event  E  as  "x  la  leas  than  or  equal  to  X"  where  X  la  a  fixed 
real  number.  Let  the  nuabar  of  trials,  out  of  the  first  H,  in 
which  E  is  observed,  be  x^,.  If  the  relative  frequency  n^/N  tends 
to  a  limit,  whatever  the  value  of  &,  that  is 


Pr  tx  <  X]  «  if* 


“h 

ST 


(1) 


exists  for  every  real  number  X,  then  x  is  called  a  random  variable. 

Of  the  multiplicity  of  random  variable  classes  two  are  presently 
of  considerable  use.  These  are  discrete  and  continuous  random 
variables.  Our  interest  herein  will  be  confined  to  the  continuous 
random  variable  which  we  conceive  as  having  the  possibility  of  taking 
any  value  over  some  interval.  More  precisely  we  want  to  be  able 
to  evaluate  the  probability 


Pr  [X, 


x  <,  X2] 


F(X) 


-  F(x) 


'x  *  X, 


(2) 


md  m  say  x  is  a  continuous  random  variable  if  F(x)  is  continuous. 
Note  that  since  F(x)  is  continuous  than  lim  Pr  [X.  <  x  <_  X-]  *4  ., 


that  is  continuous  random  variables  have  the  "odd"  property  that  the 
probability  of  taking  any  one  specified  value  is  sero. 

The  function  F(n)  is  called  a  cumulative  distribution  function. 
These  functions  have  the  following  common  properties  (s)  o  <_  F  (x)  <_lj 


(b)  F  is  a  nondecreaaing  function  of  xj  (c)  (almost  all) 

Jig  F(x)  =  1,  l^F'  )  *  0. 

If  the  function  F(x)  Is  differentiable  we  define  the  prebabiH 
density  function  p(x)  by 

p(x)  =-^  [F(x)3  (3 

and  note  that  p(x)t  like  F(x)f  represents  a  property  of  the  random 
variable  x.  Indeed  without  knowledge  of  p(x)  (or  F(x))  the  random 
variable  is  computationally  essentially  useless. 

From  the  definition  of  p(x)  it  follows  that 

X9 

Pr  [X^  <  x  <_  X^j  »  J*  p(x)  dx  (*♦ 


where  X^  < 

If  x  and  y  are  *~wo  random  variables  the  probabilities 
Pr  [x  <_  X]  and  Pr  Cy  *3  each  exist.  If  Pr  [x  X,  y  y]  exists 
for  all  X  and  all  V,  the  joint  cumulative  distribution  function  of 


x  and  y  is  defined  by  the  <.  ^oation 


F(x,y)  ix!xtyry,Prt*ilt,yiV]. 

If  both  x  and  y  are  continuous  variables  we  define  their  joint 
probability  density  function  p  (xfy)  by 


f p(x,y)  dx  dy  ■  Pr  C(x,y)  in  Rj. 


Here  R  is  any  simply  connected  region  of  the  x,  y  plane  end 


pU.y)  ■  fbjSOi 
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(b)  Basic  Properties 

Two  random  variable  x  and  y  are  mutually  Independent  if  the 
distribution  of  values  of  x  is  unaffected  by  the  value  of  y,  and 
vied  versa.  As  a  consequence  of  this  one  finds  that 

F(x„y)  =  F(x)  F(y) 5  p(x,y)  *  p(x)  p(y)  (8) 

for  independent  random  variables  x  and  y.  This  "product11  property 
greatly  simplifies  the  computational  details  for  independent 
processes. 

One  of  the  principal  features  of  a  random  process  is  its 
non-repeatability.  Therefore,  for  example,  a  time  history  of  the 
response  of  a  vehicle  to  a  random  excitation  ia  not  very  meaningful. 

Of  more  interest  is  the  probability  of  the  severity  of  the  structural 
response  to  a  given  input  occurring  throughout  the  ti’se  interval  of 
interest.  For  a  random  process  we  seek  the  mean  or  expected  value 
of  the  response.  If  x  is  a  continuous  random  variable  with  probability 
density  fun^ion  p(x)  the  expected  value  is  defined  by  the  equation 

E(x)  »  x  p(x)  dx  (9) 

where  A  and  B  define  the  range  of  values  tnet  x  can  assume.  This 
range  is  often  (—,  •)  or  (0,  -).  It  is  clear  that  £(x)  is  not  a 
function  of  xi 

If  f(x)  is  integrahle  the  expected  value  of  f(x)  is  defined  as 

ECf(x)]  *  /®  f(x)  p(x)  dx.  (10) 


4 


Sene  elementary  but  useful  properties  are  immediately  obvious.  Thus 


a)  if  a  is  constant  ECaf(x)3  3  aE[f(x)j 

and  E[f(x)  *  a]  *  E[f(x)3  ♦  a 

n  n 

&)  E[E  f,  (x)]  «  I  E[f.(x)3  (11) 

j=l  2  jsl  1 


In  particular  we  employ  the  results 


B  B 

p  (x) »*  /  y  p  (x,y)  dy,  p  (y)  =  /  *  p  (x,y>  dx. 
A  A 

y  « 


(12) 


Some  less  obvious  properties  include  the  following: 

c)  i.r  the  x^,  i  =•  lt  2,...#  n  are  mutually  independent 

then  EfS  x^,3  *  8  E  (xi). 
i=l  i«l 

d)  E  Q(x-A)r3  is  the  rth  moment  of  x  about  A.  If  A  =  0, 

E  (xr)  is  the  rth  moment  about  zero.  If  A  =  E(x)  then 

E  [(x-A)r3  -  }*  (x)  is  the  »th  central  moment.  In  terms 

r 

of  this  notation  and  E(x)  we  have 

u2  *•  var  (x)  =  o2  (x)  ■  E[(x  -  E(x))23 

3  *r  x23  -  [E(x)32  (13) 

«)  measures  the  "variability  of  x". 

p3  #3/1*2  measures  the  skewness  of  the  distribution  of  x. 

o 

measures  the  kurtosls  of  ...e  distribution 
Additional  definitions  will  be  given  where  needed  in  the  discussion. 
Details  of  analysis  and  background  are  available  in  Biller  [23  or  Parzen 
[33. 
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3. 

The  governing  equations  of  a  linear  lumped  parameter  physical 
system  are 

MX  +  CX  +  tX  3  F(t)  (14) 

where 

h  =  [M .  J  is  the  matrix  of  masses  and  inertias, 

ji 

C  3  [C . .]  is  the  matrix  of  dissipation, 

j1 

X  =  [Kjj]  is  the  stiffness  matrix, 

X  3  CXAJ  is  the  d is plan am ant  column  vector, 

F  =  [F_j]  is  the  column  vector  of  forcing  "unctions,  and 

i,  j  3  *,  £-t  •  »  . ,  n. 

Ad  is  well  known  these  equations  can  be  decoupled  with  the 
result  that 

X  3  *x  (15) 

where 

• 1  [V  15  the  matrix  of  normal  modes  (each  mode  is  a  column), 

and  x  =  [x^]  is  the  column  vector  of  normal  coordinates.  Upon  introducing 

T 

Eq.  (15)  into  Eq.  (14)  and  premultiplying  by  ♦  w  have 

^Mex  +  ♦TCdx  +  ♦^K*x  *  d^F  (16) 

IT  X 

where  ♦  M*  and  ♦  Ke  are  diagonal  and  ♦  CS  is  diagovl  if  the  damping  forces 
are  proportional  to  either  stiffness  or  mass.  In  such  cases  the  equetion 
of  motion,  Eq.  14,  can  be  written  in  the  uncoupled  form 

m5J  t  cx  +  kx  3  f  (17) 

where  m,  c  end  k  are  diagonal  matrices  and  f  »  ♦  F. 

The  response  of  an  n  degree  of  freedom  system  becomes  the  problem 
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of  finding  the  solutions  of  the  n  uncoupled  linear  differential 
equations.  If  the  individual  responses  are  desired  it  is  debatable 
whether  the  computation  should  be  carried  out  in  thi*’  form  or  in  the 
original  fora.  Often  knowledge  of  the  nortal  modes  is  highly  desirable, 
in  which  case  the  computation  as  outlined  here  is  required. 

4.  Response  of  Single  Normal  Hode  to  Deterministic  Fencing 
The  equation  for  the  jth  normal  mode  x^  is 

x.  +  2a.  w  .  x.  x.  *  f .  (t)/m.  (18, 

3  3  nj  1  nj  j  J  3 


where 


203  =  Y  (kj“j)1/2’  “nj  3  YY 


Next  a  (transfer )  function  relating  the  modal  forcing  function 
to  the  modal  displacement  x^  is  sought  in  the  fora 


If  f^  (t)  «  e^  then 


r*nj  '  “  +  2iaj“"nj] 


iwt  ,  ,  iut 

*  Hj  (»»)•  . 


The  utility  of  («)  can  be  generalized  to  any  forcing  function  that 
can  be  obtained  as  the  superposition  of  a  number  of  harmonic  forcing 
functions.  Thus  if  («)  i*  one  of  the  components  cf  the  Fourier  Series 
representation  of  f  (t)  defined  by 


Ij  (-)  fj  (t)  e”^Mtdt 

Am  the  response  toF^  (*)  is 
5^(«)  =  H^(mi  f^mi.end  the  total  response  is 
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xi (t>  ■  57  XT  ,u,t  d» 

=  J2Z  j?,  (•>  d»/_; .  it)  .*1"t  d,  (2j) 

The  function  H  (w)  is  also  the  vehicle  used  to  find  the  response  to 
a  random  excitation.  Those  relationships  are  mathematically  easier  to 
perform  ky  a  superposition  of  elemental  solutions  in  the  time  domain 
rather  than  in  the  frequency  domain.  The  relation  between  transfer 
functions  in  the  time  and  frequency  domains  is  obtained  as  follows: 

Let  5  (t-t)  be  the  unit  impulse  applied  at  t  =  t  and  h^  (t-*)  be  the 
impulse  response  function. Then  the  total'  response  to  a  continuous  farcing 
function  f  (t)  is 

xj(t)  =  (t)  hj  (t-x)  di  (23) 


Also,  we  know  from  the  linear  theory  that  the  impulse  response 
function  h^  and  («)  ar®  related  through 

hi  ^  =  §7  S_m  Hj 

(24) 

H.  (*)  =  /**  h,  (t)  e'^  dt 

1  J  —  j 


that 

This 

time 


is  they  are  Fourier  and  inverse  Fourier  transforms  of  one 
is  the  dasired  relationship  between  transfer  functions  in 
and  frequency  domains,  assuming  these  transforms  exist. 
Some  useful  transfer  functions  are  listed  here: 

(a)  Force  input,  displacement  response 

-1 

■i 

H  (M)  =  — j - - 

3  [w ,  -  ♦  2®or.  w.w] 

nj  j  nj 


another. 

the 


(25) 
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(b)  Force  input,  acceleration  response 


-1  2 

“j  (»> 5  ^ — -r“~ 


TV,  -  w  +  2ia,w  .«] 

«3  -  j  nj 


(26) 


(c)  Acceleration  input,  displacement  response 


H  («)  *  — s — — 3^ - — 

!•».-«♦  2ia.w  .w] 


nj 


j  nj 


(2?) 


(d)  Acceleration  input,  acceleration  response 

2 

H.  (»)  *  a 

3  L«.  -  «  t  2ia.« 

nj  j  nj 


(28) 


5.  Response  of  a  Single  Horaal  Mode  to  Random  Forcing 

Since  a  random  vibration  is  cot  repeatable  a  time  history  of  the 
response  of  a  structure  to  a  random  excitation  is  not  very  useful.  Of 
wore  interest  is  the  probability  of  the  severity  of  the  structural 
response  to  a  given  input  during  a  time  interval  of  interest.  Knowing 
the  relationship  between  input  and  output,  Ll(w)  or  h(t) ,  we  can  calculate 
the  statistical  properties  of  the  response  if  we  know  those  of  the  input. 

For  a  random  process  we  seek  the  sxpected  value  of  the  responao 
assuming  we  know  the  expected  value  of  the  .input.  Thus  we  seeV,  from 
Eq.  (23),  with  6  »  t  -  T 


CO]  =  E£  fZ  fj  <t-«)  hj  (0)  <163 
a  fZ  E[fj  (t-e)3  hj  (6)  d6 


(29) 


where  f ^  i»  the  only  random  variable. 
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W®  now  define  several  terns  which  are  useful  in  describing  rend  ass 
processes,  A  stationery  process  x(t)  is  one  having  statistics  which 
do  not  change  with  time.  If  its  probability  density  and  all  higher 
order  densities  are  inaependent  of  ti£%  we  call  the  process  strictly 
stationary.  An  ergodic  process  is  a  randon  process  for  which  time 


glxT 


1  /.T 


g[x(t)]  dt 


and  ensemble  averages  (another  asst*  for  expected  values)  are  equal. 
All  ergocic  processes  are  stationary  but  not  conversely. 

If  our  randan  process  is  stationary  then  it  follows  fro©  Eq.  (29) 


(t)3  =  EC fj  (t)l  f  ^  (8)  d8.  (3] 

By  setting  w«0  in  Eq.  (24)  the  value  of  the  integral  in  Eq,  (31)  can  be 
evaluated  to  obtain 


EtXj  (t)3 


Hj  (0)  ECfj  (t)3 


This  becones 


a*.  <*»  *T7 


ECfj] 


for  displacement  response  to  force  input »  from  Eq.  (?5). 

The  average  value  of  the  product  of  a  function  of  time  with  the  same 
function  displaced  t  sec  is  called  the  autocorrelation  function 
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.<s 


(34) 


R(t)  =  j  ^ItTrT  *  lia  fT  x(t)  x(t+t)  dt  ■■' 

T*-  Jo 

to  ^stlsgulsh  it  from  the  oroascorrelatien  function  g^(t)  g„  (t+t), 

R(r)  is  particularity  important  because  it  forms  the  link  with  the 
frequency  -  component  methods  of  description. 

For  an  ergodic  process  we  can  also  write  R(t)  as  an  ensemble 
average  (expected  value)  called  the  covariance  function 

R(r)  *  E[x(t)  x(t+r)3.  (35) 

2 

For  an  ergodic  process  R(t)  has  the  following  properties s  (a)  R(0)  «  x  \ 
(b)  R(-t)  a  R(t),  that  is  R  is  an  even  function*  (c)  R(0)  af  j£(.r)j  for 
all  t. 

The  relationship  between  R(t)  and  the  frequency  -  component  descrip¬ 
tion  of  a  random  function  (the  spectral  density)  is  developed  in  the 
sane  manner  as  that  relating  h^(t)  and  (m)  in  Section  4.  A  sufficient 
condition  for  the  existence  of  the  Fourier  transform  of  R(t)  is 

/_£  |r(t  )  j  dt  <  « 

Tbs  Fourier  transform  Si(w)  of  the  autocorrelation  function  of  an  ergodic 
process  is  called  the  spectral  density  defined  via 

su>  *  fl  *~Ux  dT  (36) 

(compere  with  Eq.  (24)) 
and  the  inversion  of  Eq.  (36)  is 

R(t)  =  J' -  S(w)  e^*T  <iui  (37) 

Worn  Eq.  (37)  it  is  clear  that,  the  mear.  square  value  of  the  random  variable 


a 


is  related  to  the  spectral  density  by 

7x(t)F  *  R(0)  x  J"  S(w)  d«.  (38) 

A  nuabar  of  properties  of  S(w)  are  deducifcie  free  its  definition: 

(a)  S(w)  is  even;  (b)  S(us)  is  reel  and  positive;  (c)  the  average  power 
dissipated  in  a  one  ohsi  resistor  by  those  frequency  components  of  a 
voltage  x(t)  lying  in  a  bend  between  w  and  ti  t  dm  ia  2S(w)  d» 

(the  units  of  S («)  are  power/ eycle/seeond). 

If  S^(w)  is  the  spectral  density  of  the  random  forcing  function 
fj(t)  of  Eq.  (17)  then  it  follows  (see  Aseitine  [43)  that 

$£  (w)  s  | 2  S.(w)  (39) 

4  J  1 

for  an  ergodic  pro  u«8  s^  it  follows  from  the  previous  remarks 
that  the  mean  square  response  is 

E^2]  »  [Xj(t)32  ■  R(0) 

1  *fl  V'  *• 

f.  ^7  J~  jd^(hi) | 2  S^U)  da. 

The  values  E[x^],  Eq,  (32),  and  the  root 
(*j>ts*  *  define  the  average 

The  variation  of  the  response  from  these 
the  variance,  Eq.  (13), 

02(xj)  -  E[Xj2]  -  {E[Xj 3>2, 

whose  square  root  is  the  standard  deviation. 


(u0) 

mean  square  response 
properties  of  response . 

average  values  is  given  by 
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6,  Probability  Distributions 

The  cosylete  dea-ription  of  a  .  .ndo«  process  requires  the 
selection  of  a  probebili  y  distribution.  This  is  usually  don®  by 
preserving  the  probability  density  function  p(x).  For  the  vibration 
enviroruesnt  of  vehicles  the  actual  distribution  could  be  determined 
by  extensive  processing  of  recorded  operational  data  coupled  with 
investigations  of  the  sources  of  excitation.  Investigations  of  these 
types  have  not  yet  led  to  definite  general  conclusions  on  the  nature 
of  the  probability  distributions.  For  the  purpose  of  most  structure 
and  vehicle  vibration  tests  and  responses  to  rand os  vibration  the 
Gaussian  (Horaal)  distribution  is  aest  \io«Bonly  chosen.  Its  density 
function  is 

p(x)  *  1  exnr-(x  -  !a)2/2o2]  0*1) 

err  2s 

where  «  *  Etx]  and  p2  is  the  variance. 

A  distribution  which  has  been  extensively  used  in  the  past  two 
decades  is  the  Weibull  distribution  (sse  Johnson  and  Leone  [5].  p.112) 

p(x>  *  £  (£j~  )°'1  expC-(“g— )C]  <«) 

for  x>C|  b>0(  c>0.  The  ounulative  distribution  is 

F(x)  *  1  -  exp[  -  t1*3) 

Mote  that  if  a*0,  cel  the  Weibull  distribution  includes  the  exponential 
distribution 

p(x)  «  8  expC  -  8x]f  8  a  1/b.  (*»*♦) 
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Th«  genera].  applicability  of  the  Gaussian  distribution  follows  from 
tha  Ceatral  Limit  Theorem,  Ons  of  many  forms  of  that  deHoivrc-Laplace 

h  mm  ii  rnn»e<  ■— m  mca— — m> 

result  is  as  follows:  "Let  the  random  variable  x  be  distributed  with 

mean  p  and  variance  a2  (but  with  density  function  unknown).  Then  the 

distribution  of  the  sample  mean  x  is  closely  approximated  by  the  Gaussian 

o 

distribution  with  mean  p  and  variance  a  /r.  when  n  is  large."  (Feller 
[2])  Thu*  a  process  will  be  approximately  Gaussian  whenever  the  process 
results  from  tha  superposition  of  a  large  number  of  sub-processes  in 
which  no  single  sub-process  inates  all  others. 

If  the  excitation  random  process  f£t)  has  a  Gaussian  dir+ribution 
then  the  response  x^(t)  is  also  Gaussian.  With  the  knowledge  of 
Gaussian  response  and  the  resulting  statistics  for  x^(t)  known  the 
probability  that  the  response  is  below  any  desired  level  can  be  calculated 
from  standard  tables.  For  example,  the  probability  that  the  absolute 
value  of  the  i  tsponsa  will  not  exceed  the  standard  deviation  a  is  .68, 
that  it  will  not  exceed  2o  is  .95  and  that  it  will  not  exceed  3o  is  .999. 


7 .  Typical  Calculations 

We  turn  now  to  the  evaluation  of  the  mean  response  E[x^],  and  the 
2 

mean  square  response  E[x^J ,  foe  several  typical  forcing  functions. 
Commonly,  this  input  is  specified  in  terms  of  its  frequency  limits, 
spectral  density,  and  mean  value. 

Mean  Value  of  Response 


The  mean  value  of  the  response  x^(t)  as  a  function  of  the  mean  value 

of  the  input  has  been  previously  established  ss  Eq.  (32).  Thus  we  have 
a )  mean  displacement  response  to  a  force  input: 

-  \1  1  r-,r  r  /  .  \  -I  /  .  \ 


ECx^t)]  > 


E[f j(t)] 


"j  "nj 
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b)  mean  displacement  response  to  an  acceleration  inputt 


E[x,(t)]  *  — ECA-(t)]  *46) 

J  «  f  J 

"nj 

c)  mean  acceleration  response  to  either  a  force  or  accelera  ion 
input: 


*  0. 


(47) 


Mean  Square  Response 

If  and  «2  are  the  bandwidth  frequency  limits  for  a  particular 
excitation  function  than  the  general  expression  for  the  mean  square 
response,  from  Dq.  (40),  becomes 


ECxJ'J  «  ^  J^2  iHjU)!2  gjU)  da 


(48) 


for  the  esse  ft^(«)  *  ■  constant  (‘t^a/cps) 


Etx^] 


■-^?r  - 

2*  m.  (« 


hJ 

ds 


,  “5  2  vi  2  2  2 

■J  "I  Un]  -  «  >  *°3  “nj  “ 


(49) 


is  the  mean  square  displacement  response  to  a  constant  force  input. 

B/3 

If  tha  excitation  ia  of  tha  form  8^(w)  *  Aoj^  * 

<_  w  <_  »2  (Ihsftbpa)  then 


EtX^] 


B/3 

•H  d« 


,  H  .  2  ,  2  .  2  » 

1  ♦  2.  .  U  -1)  .  *  «Bj 


(50) 


is  the  swan  square  acceleration  responsa  to  a  varying  force  input. 
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Similar  express  Irr.®  arc  obtainable  for  other  mean  square  response*. 
The  .integral a  can  be  evaluated  in  t«m«  of  elementary  functions. 

8.  Response  of  a  Complax  Struetur  to  Random  Excitation 

The  solutions  of  the  random  responses  of  each  nomal  node  must  now 
be  combined  to  yield  the  random  response  of  the  multi-degree  of  freedom 
system.  Modal  (f^)  and  physical  forcing  functions  (F^)  are  related 
by  the  equation  (see  Eq.  (16)) 


f3  *  k  k  ri 


where  4.^  are  the  elements  of  4  .  Thus  we  have 


(51) 


=  ll  ♦ji  ECFi3  (52) 

that  is  the  mean  value  of  the  modal  forcing  function  is  directly 
computable  in  terms  of  the  mean  values  of  the  physical  forcing 
functions.  Equation  (52)  again  follows  from  the  linearity  of  the 
expectation  operator. 

The  relationship  between  the  physical  and  aodal  spectral 
densities,  used  in  our  previous  work,  is  more  complex.  From  specific 
font  of  Eq.  (38)  the  relationship  betwesn  input  force  and  input 
spectral  density  is 


ELFi2(t)]  »  ^ 

and  in  the  model  plane 

ECfj2(t>]  .  ^ 


"5^  (w)  dw 


Sj  (*»)  d*j 


(53) 


(54) 
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aisSKsyi! 


Since 


f  2  =  §  ♦.  I  F 

j  r.,  8=1  l*  js  r  8 


then 


ECf^  -  L  *jr  ECFr2]  *  L  LV<»  Fr  F» 

Using  F.qs.  (55)  and  (55)  it  follows  from  Eq,  (40,,  for  ergodic 
that 


SjU)  *  ?, a  *jr  *js  CV“  ?8Uie8);1 

where  6  is  the  phase  angle  of  the  excitation. 

9.  From  Modal  Response  to  Physical  Response 
In  general  X  =  #x  or  in  component  farm 


1/2 


Xi  * 


?=i  ♦«  Y 


Consequun  'y  the  mean  physical  response  is 


w  *  ?.i  *ij  “V 


and  the  mean  square  response  is 


EtXi2]  *  ♦  2  E[Xj2]  +  ^  4ik  ECxjX)cJ. 

if  and  x^  are  independent  the  quantity  ECx^x^}  »  cov  (x^.x^i 
in  which  case 

X*L*'  -  j=1  ^  £t*j23. 


(55) 

(56) 

processes , 

(57) 

(58) 

(59) 

(60) 

1  *  °* 

(61) 
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This  situation  occurs  sufficiently  often  to  be  of  interest  in 


real  applications.  Finally, 

o*  (Xi)  =  ECX^]  -  {ECX^]}2.  (62) 

Since  the  relationship  between  physical  and  nodal  acceleration  is 
•*  n  ,, 

V  5=1  ♦«  h 

the  steps  outlined  above  are  repeatable  for  the  acceleration  response 
of  the  linear  system. 

Other  statistical  considerations  also  play  an  important  role  in 
this  Themis  Project  and  therefore  will  be  the  subject  of  additional 
summary  reports.  These  include  fatigue  damage  criteria  for  structures 
under  random  excitation,  the  response  of  oontinuoua  (rather  than  lumped) 
atructur#to  random  excitation  and  the  optimisation  of  lumped  parameter 
and  continuous  dynamic  systems  under  random  excitation. 
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